Superconducting phase diagram of the filled skuterrudite PrOs 4 Sbi 2 
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We present new measurements of the specific heat of the heavy fermion superconductor PrOs4Sbi2, 
on a sample which exhibits two sharp distinct anomalies at T cl = 1.89K and T C 2 = 1.72K. They 
are used to draw a precise magnetic field-temperature superconducting phase diagram of PrOs4Sbi2 
down to 350 mK. We discuss the superconducting phase diagram of PrOs4Sbi2 and its possible 
relation with an unconventional superconducting order parameter. We give a detailed analysis of 
H C 2(T), which shows paramagnetic limitation (a support for even parity pairing) and multiband 
effects. 

PACS numbers: 65.40.Ba,71.27.+a,74.25.Dw,74.25.Op,74.70.Tx 



I. INTRODUCTION 

The first Pr-based heavy fermion (HF) superconductor 
PrOs4Sbi2 (T c ~ 1.85K) has been recently discovered^. 
Evidence for its heavy fermion behavior is provided 
mainly by its superconducting properties, like the height 
of the specific heat jump at the superconducting tran- 
sition or the high value of H C 2(T) relative to Tc 1 . 
PrOs 4 Sbi2 is cubic with 2), point group symmetry, and 
has a nonmagnetic ground state, which in a single ion 
scheme can be either a T23 doublet or a Ti singlet, a ques- 
tion which remains a matter of controversy. Presently, 
most measurements in high field seem to favor a singlet 
ground state 3 *^^. In any case, whatever the degeneracy 
of this ground state, the first excited state (at around 6K) 
is low enough to allow for an induced electric quadrupo- 
lar moment on the Pr 3+ ionsi&, that could explain the 
heavy fermion properties of this system by a quadrupo- 
lar Kondo effect. Thus, while the pairing mechanism of 
usual HF superconducting compounds (U or Ce-based) 
could come from magnetic fluctuations, the supercon- 
ducting state of PrOs 4 Sb 12 could be due to quadrupolar 
fluctuations. Yet at present, this attractive hypothesis 
is backed by very few experimental facts, both as re- 
gards the evidence of a quadrupolar Kondo effect in the 
normal phase and as regards the pairing mechanism in 
the superconducting state. Even the question of the un- 
conventional nature of its superconductivity is still open. 
Indeed, several types of experiments have already probed 
the nature of this superconducting state, but with appar- 
ently contradictory results. Concerning the gap topology, 
scanning tunneling spectroscopy measurements point to 
a fully open gap, with some anisotropy on the Fermi 
surface 9 . Indication of unconventional superconductivity 
might come from the distribution of values of the resid- 
ual density of states (at zero energy) on different parts of 
the sample surface. This could be attributed to a pair- 
breaking effect of disorder. The same conclusion as re- 



gards the gap size was reached by ixSR measurements 10 
and NQR measurements 11 , although unconventional su- 
perconductivity is suggested in the latter case by the ab- 
sence of a coherence peak below T c in \/T\. 

This should be contrasted with recent penetration 
depth measurements, that would indicate point nodes of 
the gapi£, or the angular dependence of the thermal con- 
ductivity which suggests an anisotropic superconducting 
gap with a nodal structure^. This latter measurement 
also suggests multiple phases in the temperature (T) - 
field (H) plane, which could be connected to the dou- 
ble transitions observed in zero fielcU^i^*i&. Recent fiSR 
relaxation experiments^ detected a broadening of the 
internal field distribution below T c , suggesting a multi- 
component order parameter or a non unitary odd parity 
state, with a finite magnetic moment. 

In this context, our results bring new insight on the 
question of the parity of the order parameter, and draw 
a definite picture of the (H,T) phase diagram as deduced 
from specific heat measurements. With reference to the 
historical case of UPt3, we emphasize that the present 
status of sample quality may explain the discrepancies 
between the various measurements : definite claims on 
the nature of the superconducting state in PrOs4Sbi2 
are at the very best too early, the key point being the 
sample quality. 



II. EXPERIMENTAL DETAILS 

We present results on single crystals of PrOs4Sbi2 
grown by the Sb flux methodi^*^. These samples are 
aggregates of small single crystals with well developed 
cubic faces. They have a good RRR of about 40 (between 
room temperature and 2K), a superconducting transition 
(onset of C p or p) at 1.887 K and, they present a very 
sharp double superconducting transition in the specific 
heat. 
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FIG. 1: Specific heat of samples of the same batch including 

Q 

sample n°l as versus T at zero field measured with a 
quasi-adiabatic method. The inset is a zoom on the double 
superconducting transition of sample n°l measured with an 
ac method: T cl = 1.887 K and T c2 = 1.716 K. 
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FIG. 2: Temperature sweeps of the ac specific heat of sample 
n°l at several fields below 1.4 T. The normal state was sub- 
tracted. The arrows indicate the double transition T c i and 
T c2 . 



Two different techniques have been used for the specific 
heat measurements. The first is a quasi-adiabatic method 
with a Au/Fe-Au thermocouple controlled by a supercon- 
ducting quantum interference device (SQUID) in a 3 He 
calorimeter. It is well suited to quantitative studies in 
zero field (the addenda are precisely known), and was 
used on samples with a total mass of 8 mg. The second 
technique was ac calorimetry, used to follow the super- 
conducting transitions under magnetic field in order to 
draw a complete phase diagram of the two superconduct- 
ing transitions. This ac calorimetry uses a strain gauge 
heater (PtW alloy), a sensitive SiP thermometer (silicon 
doped with phosphorus close to the critical concentra- 
tion of the metal- insulator transition), and a long gold 
wire (25/im diameter) as a heat leak. For the ac method, 
we choose a frequency of 0.04 Hz and a large integration 
time of 350 s. The SiP thermometer is measured with a 
four lead resistance bridge at 500 Hz, whose analog out- 
put is sent to the lock-in detection. The heating power 
was chosen so that the SiP temperature oscillations re- 
main smaller than 6 mK in order to avoid broadening of 
the transitions. Thermometry under field was controlled 
by thermometers located in the (zero field) compensated 
region of the magnet. 



III. SPECIFIC HEAT RESULTS 

Figure^displays the specific heat C p (0, T) /T of 3 sam- 
ples of the same batch measured together, of total mass 
8 mg. The inset of fig. Q]is the ac specific heat C p (0,T) 



of one of these three samples (thereafter called sample 
n° 1). As it has been previously observed^, PrOs4Sbi2 
shows a Schottky anomaly with a maximum in C p /T 
near T — 2.2K. Absolute values of C p /T are the high- 
est ever reported: at T = 1.7K, C p /T = 3.65J/ K 2 .mol 
and at T = 2K, C p /T=2.9 J/mol.K 2 . Sample n°l has a 
well defined double transition: to our knowledge, it is the 
sharpest ever reported in the literature, although we are 
aware of similar (yet unpublished) results by Y. Aok: 20 
on samples grown in the same group. The width of the 
two transitions was estimated to be 16mK and 58mK, 
with respectively T c2 = 1.716K and T cl = 1.887K ( with 
the junction criterion). 

Shown in fig. and fig. El are the ac specific heat 
measurements C P (H,T) at, respectively, constant mag- 
netic field (T c \ > 1.16 K) and constant temperature 
{T c i < 1-15 K) of the same sample n°\. The normal 
state specific heat, or an arbitrary line between the two 
transitions, has been subtracted for the temperature or 
field sweeps respectively. Even under field, the transi- 
tions remain very sharp, so that we were able to detect 
them down to 350 mK and to draw a precise phase dia- 
gram (Fig. 2J. The width of the two transitions does not 
exceed AH c2 = 80 mT and AH' = 50 mT at 500 mK. 

The phase diagram has the same features as reported 
by Tayama et al£ from magnetization measurements. 
The advantage of specific heat is to give an unambigu- 
ous signature of a bulk phase transition, that cannot be 
confused with other physical phenomena like peak-effect. 
The two transition lines remain almost parallel and we 
will see that they can be deduced from each other simply 
by scaling T c . 
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FIG. 3: Field sweeps of the ac specific heat of sample n°l at 
several temperatures. An arbitrary line between the two tran- 
sitions was subtracted. We follow the two transitions (_ff C 2 
and H ) down to 350 mK. 
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IV. SAMPLE CHARACTERIZATION 

Three pieces of sample n°l have been used for further 
characterizations, called la, lb and lc. As well as the 
specific heat of sample n°la (2 mg), we have measured 
the resistivity of samples n°lb (0.2 mg) and n°lc, and the 
ac susceptibility and dc magnetization of samples n°la 
and n°lb. 

Concerning the specific heat, the high absolute value 
of C p as well as the large height of the two supercon- 
ducting jumps (A(C P /T) = 350 mJ/mol.K 2 at T cl and 
A(C P /T)=300 mJ/mol.K 2 at T c2 ) must be linked to the 
high quality of these samples (absence of Sb-flux and/or 
good stoichiometry) . Moreover, the heights of the two 
steps of sample n°la are quantitatively similar to those of 
the entire batch (7.5 mg) as Vollmer has already pointed 
ou04 

Like in previous work 1 ^^, we have noticed that the 
resistivity at 300K is very sample dependent ( from 
200 to 900 /iSl.cm). On all samples, the value of p at 
300 K seems to scale with the slope at high temperature 
(T > 200K), i.e. the phononic part of the resistivity, as if 
the discrepancies were due to an error on the geometric 
factor. This error could be explained by the presence of 
microcracks in the samples. We have taken this problem 
into account by normalizing all data to the slope dp/dT 
at high temperature (T > 200K) of sample n°lc, chosen 
arbitrarily. For samples n°lb and n°lc respectively (Fig. 

, the RRR (ratio between 300K and 2K values) are 44 
and 38, the onset T c are 1.899 K and 1.893 K (match- 
ing the critical temperature obtained by specific heat), 
and the temperatures of vanishing resistance (Tr = o) are 
1.815 K and 1.727 K. Tr = o of sample n°lc is equal to 
T C 2 and this remains true under magnetic field. So, in 
sample n°lc, the resistive superconducting transition is 
not complete between T c \ and T c i. 

FigurcElshows the superconducting transition for sam- 
ples n°la and n°lb by ac-susceptibility (H ac = 0.287Oe), 
corrected for the demagnetization field. The onset tem- 
perature is the same for the two samples (1.88K). The 
transition is complete only at around 1.7K and two tran- 
sitions are visible. The field cooled dc magnetization of 
samples n°la and n°lb (Hd c = lOe), shown in fig[7| gives 
a Meissner effect of, respectively, 44% and 55%, indicat- 
ing (like specific heat) that the superconductivity is bulk. 
The two transitions are also visible. 



V. DISCUSSION 



FIG. 4: H — T superconducting phase diagram of PrOs4Sbi2 
determined by specific heat measurements on sample n°l. 
The field dependence of T c i and Tci are completely similar. 
The lines are fits by a two-band model of the upper critical 
field (Section IWEH . Only T c has been changed from H C 2 to 
H . 



A. Double superconducting transition 

Let us first discuss the nature (intrinsic or not) of 
the double transition observed in our specific heat mea- 
surements. The remarkable fact, compared to previous 
reportsi^i^, is the progress on the sharpness of both tran- 
sitions. If for previous reports, a simple distribution of 
T c due to a distribution of strain in the sample could 



1000 

E 

o 

d 



sample nib 
sample nlc 




1.7 1.75 1.8 1.85 T(K) 



50 100 150 T(K) 300 



FIG. 5: Resistivity of samples n°lb and n°lc normalized to 
the slope at high temperature of sample n°lc. The inset is 
a zoom on the superconducting transition. The resistivity of 
sample n°lc is zero only at T C 2. 
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FIG. 7: DC magnetization of sample n°la at Hoc = lOe 
with zero field cooled (ZFC) and field cooled (FC) sweeps. 
The Meissner effect is 44% for this sample and 55% for sample 
n°lb (not shown). The superconductivity is bulk. 
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FIG. 6: Real part \' of the AC susceptibility of samples n°la 
and n°lb measured with an AC magnetic field of 0.287Oe 
at 2.11Hz. Like in the results of resistivity measurement, the 
superconducting transition is not complete at T c \ and the two 
transitions are visible. 



have explained the transition width, it is not the case 
anymore for the sharp features observed on these new 
samples. Also the explanation recently proposed^ that 
the lower transition at T c2 would be induced by Joseph- 
son coupling of one sheet of the Fermi surface to another 
with transition temperature T c \ is excluded, owing to the 



sharpness of both transitions and particularly of that at 
T C 2'. see the broadening calculated by the authors^ for 
the Josephson induced transition. Nevertheless, some of 
our results still cast serious doubts on the intrinsic na- 
ture of the double transition. Indeed, the susceptibility 
also shows a " double" transition, and resistivity becomes 
zero only at T C 2- If the first transition at T c \ was bulk- 
homogeneous, the resistivity should immediately sink to 
zero below T c i, and the susceptibility (x) should also 
show perfect diamagnetism far before T c i : the sample 
diameter is at least 1000 times larger than the penetra- 
tion depth (A) so that the temperature dependence of 
A cannot explain such a transition width of x (contrary 
to the statements of E.E.M. ( 'hi;: - ). So both resistiv- 
ity and susceptibility are indicative of remaining sample 
inhomogeneities. 

Nevertheless, in our opinion, even the comparison 
of the various characterizations of the superconducting 
transition by resistivity, susceptibility and specific heat 
on the same sample does not allow for a definite conclu- 
sion. 

Two historical cases are worth remembering. URu2Si2 
showed a double transition in the specific heat in some 
samples, with inhomogeneous features detected in the 
susceptibility 23 . In that case, the authors of reference^ 
could clearly show that it was not intrinsic (maybe aris- 
ing from internal strain ?) because different macroscopic 
parts of the same sample showed one or the other tran- 
sition. In PrOs 4 Sbi2, the specific heat results are repro- 
ducible among various samples of the same batch (see 
samples n°l and n°la of the present work), and such an 
easy test does not work. 
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The second case is of course UPt3: it is now well estab- 
lished that the two transitions observed in zero field are 
intrinsic and correspond to order parameters of different 
symmetries. But the first results on samples that were 
not homogeneous enough showed exactly the same be- 
haviour as the present one in PrOs4Sbi2: two features in 
the susceptibility and a very broad (covering both tran- 
sitions) resistive transition 24 . It was not until the sample 
quality improved significantly that resistivity and suscep- 
tibility transitions matched the higher one 2 ^. The puz- 
zling result for PrOs4Sb 12 , compared to UPt 3 , is that 
despite the sharpness of the specific heat transitions, re- 
sistivity and susceptibility reveal inhomogeneities, which 
means that this new compound probably has unusual 
metallurgical specificities. 

Continuing the parallel with UPt3, basic measure- 
ments rapidly supported the intrinsic origin of the two 
transitions: they were probing the respective field and 
pressure dependence of both transitions. Indeed, a com- 
plete (H-T) phase diagram was rapidly drawn, showing 
that in UPt^, like in U\- x Th x Be\^, the two transi- 
tions observed in zero field eventually merged under mag- 
netic field, due to a substantial difference in dT c /dH. It 
is even more true for the pressure dependence of T c i.2, as 
the thermal dilation has jumps of opposite sign at the two 
transitions, indicating opposite variations of T c \ 2 under 
pressure (Ehrenfest relations) 28 . So in this compound, 
the two transitions could be rapidly associated with a 
change of the symmetry of the order parameter indicating 
the unconventional nature of the superconducting state. 

We are not so lucky in the case of PrOs4Sbi2: indeed, 
the field dependence of T c i seems completely similar to 
that of T c i (fig. 0J, a claim that will be made quantita- 
tive below. It is the same situation as in UPa^Sia 2 -^, and 
nothing in favor of an intrinsic nature of the double tran- 
sition can be deduced from this phase diagram. Another 
phase diagram has already been established by transport 
measurements, with a line H*(T) separating regions of 
twofold and fourfold symmetry in the angular depen- 
dence of thermal conductivity under magnetic field 13 . 
It may seem likely^ 3 ^ that this line would merge with 
the double transition in zero field. From the line H (T) 
drawn from our specific heat measurements (T C 2(H), fig. 
U), we can onclude that this is not the case : H (T) does 
not match the line H*(T) drawn in reference^, unless 
there is an unlikely strong sample dependence of these 
lines. 

Comparison of the pressure dependence of T c \ and T C 2 
seems more promising: contrary to case of UPtj^ the 
jump of the thermal expansion at the two superconduct- 
ing temperatures does not change signAS, but from the 
relative magnitude of these jumps and our specific heat 
peaks, we get a value dTci/dp « — 0.2K/GPa, and twice 
as much for dTc2/dp. This supports a different origin 
for both transitions. The weak point is that the ther- 
mal expansion measurements were done on two samples 
mounted on top of each other, that were early samples 
with rather broad specific heat transitions. Thus, the 



question of the intrinsic nature of the double supercon- 
ducting transition remains open. 



B. Upper critical field 

Another quantity which has not been thoroughly dis- 
cussed is the upper critical field H C 2 (T) . In heavy fermion 
superconductors, H C 2(T) has always proved to be an in- 
teresting quantity, mainly due to the large mass enhance- 
ment of the quasiparticles. Indeed, the usual orbital lim- 
itation is very high in these systems, due to the low 
Fermi velocity, so that the authors of referencei could, 
from their H C 2 (T) data, confirm the implication of heavy 
quasiparticles in the Cooper pairs (revealed also by the 
specific heat jump at T c ). They also gave an estimate of 
the heaviest mass: ~ 50 mo, where mo is the free electron 
mass. 

But, as the orbital limitation is very high, H C 2(T) is 
also controlled by the paramagnetic effect. A quantita- 
tive fit of H C 2(T) easily gives the amount of both limi- 
tations, except that on this system, H C 2(T) has an ex- 
tra feature: a small initial positive curvature close to 
T c . This feature has been systematically found, what- 
ever the samples and the techniques used to determine 
H c2 (T) (p, x and Cp) 1 ' 13 ! 14 ' 21 ■ Our data, obtained by p 
on sample n°lc and C p on sample n°l, matches all other 
published results. So we can now consider this curvature 
of H C 2(T) as intrinsic, and not due to some artifact of 
transport measurements, or coming from inhomogeneity 
in the sample. Such an intrinsic positive curvature also 
appears in MgBj^ii or in borocarbides like YM2B2C and 
LuNiaBaC 3 !. 



1. Physical inputs 

We propose an explanation based on different gap am- 
plitudes for the different sheets of the Fermi surface of 
PrOs4Sbi2, which is made quantitative through a "two- 
band" model 3 ^ . Microscopically, STM spectroscopy also 
reveals an anisotropic gap, with zero density of states at 
low energy (fully open gap) but a large smearing of the 
spectra™. Recent microwave spectroscopy measurements 
also discuss a two band model^ 2 ., but in order to explain 
the double transition. We insist that our model has noth- 
ing to do with the double transition, which clearly in- 
volves heavy quasiparticles both at T c \ and at T C 2 (see the 
size of the two specific heat jumps): our aim is a quan- 
titative understanding of H C 2(T), based on the normal 
state properties of PrOs4Sbi2, as in MgB2 or borocar- 
bides where no double transition has ever been observed. 
The physical input of a multi-band model for H C 2(T) is 
to introduce different Fermi velocities and different in- 
ter and intra band couplings. As a result, T c is always 
larger than for any of the individual bands 33 . The slope 
of H C 2 at T c is larger for slower Fermi velocity (heavier) 
bands. Positive curvature of H C 2(T) is easily obtained 
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if the strongest coupling is in the heaviest bands (large 
H c2 ), with a slight T c increase due to the inter band cou- 
pling to the lightest bands (small initial slope)2i. 



2. The two-band model: 

There are at least three sheets for the Fermi surface 
of PrOs4Sbi2, but in the absence of a precise knowl- 
edge of the pairing interaction, a full model would be 
unrealistic, having an irrelevant number of free parame- 
ters. A two band model is enough to capture the physics 
of anisotropic pairing, although only the correspondence 
with band calculations becomes looser. In our model, 
band 2 would correspond to the lightest (/3) band de- 
tected by de-Haas van Alphen measurements, and band 
1 would be a heavy band having most of the density of 
states. Indeed, the de Haas-van Alphen experiments on 



PrOs 4 Sb] 



, 34.35 



reveal the presence of light quasiparti- 



cles (band (3) and heavier particles (band 7). The heavi- 
est quasiparticles are at present only seen by thermody- 
namic measurements (C p or H C 2). Anisotropic coupling 
between the quasiparticles is considered in the framework 
of an Eliashberg strong coupling two-band model^ in the 
clean limit, with an Einstein phonon spectrum (charac- 
teristic "Debye" frequency O). Let us point out that the 
results do not depend on (and a fortiori do not probe) 
the pairing mechanism, which is likely to be much more 
exotic than the usual electron-phonon mechanism. Com- 
pared to a single band calculation, there is now a matrix 
of strong coupling parameters Ajj describing the diffu- 
sion of electrons of band i to band j by the excitations 
responsible for the pairing. 

What matters for H C 2{T) is the relative weight of the 
Xij , not their absolute value: we consider T c or the renor- 
malized Fermi velocities as experimental inputs. Ajj de- 
pends both on the interaction matrix elements between 
bands i and j, and on the final density of states of 
band In MgB2, it is claimed that electron-phonon 
coupling is largest within the a bands. Here, knowing 
nothing about the pairing mechanism, we assume con- 
stant inter and intra band coupling, so that the rela- 
tive weight of the \ t j is only governed by the density 
of states of band j. This density of states is itself pro- 
portional to the contribution of that band to the specific 
heat Sommerfeld coefficient: 500mJ/K 2 .moli for band 
1, 20mJ/K 2 .mol for band 2 34 i 35 . The Fermi velocity of 
band 1 (not observed by de Haas-van Alphen measure- 
ments) is the main adjustable parameter of the fit : we 
find Dpi = 0.0153 10 6 m/s _1 , in agreement witbA where 
the Fermi velocity has been inferred from the slope of 
H C 2(T) at T c ignoring the initial positive curvature. All 
other coefficients are either arbitrary (Ai i=l) (in agree- 
ment with the strong coupling superconductivity con- 
cluded inii), conventional values (gyromagnetic ratio for 
the paramagnetic limitation g = 2, Coulomb repulsion 



parameter /x* ■ = OAS, 
(T c = 1.887 'K=> fl 



, or fixed by experimental data 
21. 7K, v F2 = 0.116.1Q 6 m.s _1 
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FIG. 8: Open circles show the data of Hca(T) by specific heat 
measurement on PrOs4Sbi2 (sample n°l) . The lines show fits 
with a two-band model (solid line), with a single-band model 
(dashed line), without the paramagnetic limit (g=0) (dotted 
line). It shows that the increase of T c due to the coupling 
with light qp band is rapidly suppressed in weak magnetic 
fields. H C 2 is also clearly Pauli limited, supporting a singlet 
superconducting state. The parameters for the solid line fit 
are: g = 2, n*j = O.l^j, v F1 = 0.0153 10 6 m/s _1 (heavy 



band), 



VF2 



0.116.10 m.s (lightest band, from de Haas- 



van Alphen oscillations 



T c = 1.887 K= 



fi = 21. 7K. 



from the de Haas- van Alphen data on the (3 band). The 
model fits well the experimental data (cf. fig. , includ- 
ing the small positive curvature. Before discussing the 
interpretation of the fit as regards the values of the pa- 
rameters and the parity of the superconducting order pa- 
rameter, let us note that we can fit the H (T) line (fig. 0J) 
with the same parameters as for H C 2 except Q, adjusted 
to give T c — 1.716K. There is a good agreement with all 
data except at very low temperature or near T c where 
the curvature is reduced compared to H C 2(T). However, 
these deviations are weak, and this is why we claim that 
H (T) has the same behavior than H C 2(T), which does 
not help to identify the second transition with a symme- 
try change of the superconducting order parameter. 



3. Interpretation 

Shown also in fig. |H1 are the calculations of H C 2 (T) for 
a single band model: Vfi, the characteristic frequency 
f2 and Ai^ have the same values as before, but all other 
Aij coefficients have been turned down to zero, elimi- 
nating the effects of the light electron band. T c is then 
reduced (down to 1.78 K) and the positive curvature dis- 
appears. We also observe that the fit of H C 2 is basi- 
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cally unchanged above 1 T, meaning that low fields sup- 
press the superconductivity due to the light electron band 
restoring a "single band" superconducting state. This 
is the same effect observed more directly in MgB2 with 
specific heat measurements under magnetic field : the 
smaller gap rapidly vanishes, leading to a finite density 
of states at the Fermi level under magnetic fields due 
to the 7r band 3 -. This suppression of the light quasi- 
particle superconductivity would have here an effect on 
specific heat too small to be observed (contribution of 
the light quasiparticles to the specific heat of order 4% 
of the Sommerfeld coefficient, itself buried in the large 
Schottky anomaly). But it may have much larger ef- 
fects on transport. Let us note that the clean limit is a 
posteriori justified: from vpi, we find a coherence length 
£o ~ llOA, whereas from a residual resistivity po and spe- 
cific heat coefficient 7 ~ 0.5J/K 2 .mol- 2.10 3 J/K 2 .m 3 , 
3-L 

we get vfiI ~ ~ 2.10 _3 m 2 /s yielding a mean free 

PI 

path I ~ 13001 > 10£o- 

More interestingly, the quantitative fit of H c2 also al- 
lows a discussion of the parity of the order parameter 
in PrOs4Sbi2. Indeed, like other heavy fermion super- 
conductors, despite the low-T c value, H C 2 can be sen- 
sitive to the Pauli limit in case of singlet pairing. The 
fit in fig. [S] includes this paramagnetic limitation, with 
the conventional free electron value for g (g = 2). Also 
shown in fig. |HJis the calculation of H C 2(T) with the same 
parameters but for g = 0, i.e. without any paramag- 
netic limit. The strong deviations observed demonstrate 
that the paramagnetic effect controls H C 2 at low T in 
PrOs4Sbi2. Quantitatively, the paramagnetic effect de- 
pends on the coupling strength. Yet, we choose arbitrar- 
ily A = 1. A rather strong coupling regime is supported 
by NQR experiments 11 . Even for a weak coupling pic- 
ture (A = 0.6), the fit yields g = 1.55. In both cases, 
the paramagnetic limit remains important, supporting a 
singlet nature of the superconductivity, contrary to that 
had been suggested io2L2£i2£ and in agreement with the 
supposition iriSS. This result should be quite robust, 
independent of the two band model. It could be inval- 
idated if the mass renormalization mechanism was field 
dependent, which could be an explanation for the dif- 
ference in the large specific heat Sommerfeld coefficient 



obtained in low field, and the de Haas-van Alphen mea- 
surements performed at high field^. In such a case, the 
"saturation" of H c2 (T) at low temperature could arise 
from a reinforcement of the orbital limitation alone. 



VI. CONCLUSION 

To conclude, we have drawn a very precise supercon- 
ducting phase diagram of PrOs4Sbi2 down to 350 mK. 
by a specific heat measurement. We have yet no clear ev- 
idence of the unconventional nature of the superconduct- 
ing order parameter from this phase diagram. The su- 
perconducting phase diagram with the symmetry change 
of the order parameter drawn by K. Izawa et ali£ from 
transport measurements does not seem related to the 
double transition observed with specific heat measure- 
ments. Despite the high quality of the sample, we can- 
not completely exclude that there are still two parts with 
different T c in our sample, as H (T) is just scaled from 
H C 2 with respect to T c . The puzzling result is that de- 
spite sharp specific heat transitions, inhomogeneities are 
still present in the samples. This calls for caution in the 
claim of various types of nodes of the gap by different 
sophisticated techniques: the most urgent task is to un- 
derstand the problem of sample quality. Contrastingly, 
the upper critical field is very reproducible, independent 
of samples and types of measurements. It has been ana- 
lyzed with a strong coupling two-band model taking ac- 
count of the spread in the effective masses of the quasi- 
particles and of the pairing strength as suggested also by 
STM spectroscopy measurements^. The strong influence 
of the paramagnetic limit on H c2 is the first experimental 
argument for a singlet superconducting order parameter. 
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